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ABSTRACT 


An attempt to describe the atmospheric boundary layer during 
the MIZEX 84 from 11 July to 17 July was made using a sodar 
system. The computation of the temperature structure parameter CT? 
in the surface layer using in situ measurements allowed the calibra- 
tion of the system. CT? was found to have the following functional 
dependence on the backscattered signal I: 

Log(CT?)= 8.63Log(I) -25.66. 

Time-height cross-sections showed very well the structure of 
the boundary layer. An especially good representation of the 
strength of CT? in the inversion layer was achieved. 

A program was developed to obtain the inversion layer height 
Zi, where the signal reaches a maximum. The thickness of the 
inversion layer was also computed and was used with the calibration 
law to compute CT? at the inversion layer. This allowed the 
computation of the jump of potential temperature and the refractive 
index gradient in the inversion under free convection. Refractive 
trapping conditions were evident on 15 July although the generally 
small value of the heat flux did not indicate a strongly unstable 
boundary layer. Comparaison with two  radiosonde profiles shows 
good agreement in one case, but does not permit a reliable conclu- 
sion to be made concerning the method. However, the procedure 
could be extended to other areas, especially those with strong free 
convection. The method should provide a good estimate of the 


refractive condition in the atmospheric boundary layer. 


RESUME 


Un systeme sodar (Sound Detection and Ranging) a permis de 
recueillir des données concernant la couche limite atmosphérique 
entre le 11 et 17 juillet pendant l'opération MIZEX 1984. Le calcul 
du paramètre de structure de la température CT? dans la couche de 
surface utilisant des mesures in situ, permet la calibration du 
systeme. Il apparait que CT? est fonction du signal de reverbera- 
tion direct I suivant la loi: 

log(CT?) = log(I) - 25.66. 

Plusieurs coupe représentant CT? en fonction du temps et de 
l'altitude montrent la structure de la couche limite et donnent une 
bonne représentation de la valeur de CT? dans la couche 
d'inversion. 

Un programme a ete développe pour obtenir la hauteur de cette 
couche ou un signal maximum est atteint. L'epaisseur de cette 
couche est aussi calculée afin d'etre utilisée avec la loi de calibra- 
tion pour en déduire CT? a chaque niveau. Ceci permet de 
trouver le saut de temperature potentielle et le gradient d'index de 
refraction dans l'inversion en cas de convection libre. Pendant la 
journée du 15 juillet on remarque que les ondes electromagnetiques 
ont pu être "capturees" par la couche d'inversion, bien que la 
faible valeur du flux de chaleur n'ait pas ete toujours conforme aux 
conditions d'instabilite requises. La comparaison avec quelques 
profiles de radiosonde est bonne dans un cas, mais ne permet pas 
d'obtenir une conclusion certaine sur la methode. Il pourrait cepen- 
dant étre interessant de l'étendre à d'autres zones soumises à une 
plus forte convection libre, ceci donnant une estimation des condi- 


tions de reéfraction dans la couche limite atmosphérique. 
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I. INTRODUCTION. 


The atmospheric boundary layer (ABL), which is found just 
above the surface and generally extends to several hundred meters, 
can be considered to be composed of three different layers. At the 
lowest level is the surface layer, whose height is about one tenth 
of the ABL height and where turbulent exhanges with the surface 
are driven directly by vertical gradients. The mixed layer comprises 
most of the remaining ABL depth. Under unstable conditions 
mechanical and thermal mixing cause homogenous conditions. The 
third layer is the inversion layer where turbulent kinetic energy is 
small and there are strong vertical gradients. In particular, the 
gradient of potential temperature is associated with stronger fluctua- 
tions of the temperature. 

The polar region is of particular interest because it is a 
region not only of scientific interest but also of military importance. 
This area differs significantly from the midlatitude and tropical 
ABL's, where anomalous electromagnetic propagation due to mean 
vertical gradient of temperature and humidity are frequently 
observed. Although anomalous propagation is not commonly reported 
in the Arctic Ocean, the boundary layer can extend to a low 
height and have a strong inversion layer. Another difference 
concerns the difference in the stability condition, with a generally 
more stable atmosphere at high latitudes and more convective condi- 
tions in mid-latitude or in the tropics. Therefore, it is important 
to understand the vertical structure of Arctic atmospheric boundary 
layer. To date, little work have has been done on this problem. 

The Marginal Ice Zone experiment (MIZEX 84) was conducted in 
the summer of 1984 and produced a wealth of data with which to 
investigate the structure of the Arctic ABL. Several instrumented 


platforms were deployed in the East Greenland Sea during that 
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summer to measure the mesoscale processes that influence the 
advance and retreat of the ice margin. A variety of data were then 
gathered, including measurements of surface variables and  sodar 
data. The sodar (Sound Detection and Ranging) system is one of 
the sensors used in this study. It measures the strength of the 
small-scale turbulence encountered along its vertical ray path. The 
intensity of the signal is proportional to the refractive-index struc- 
ture parameter Cn’, called also the structure function parameter, 
which has an important environmental influence on the propagation 
of electromagnetic waves. Profiles of this parameter provided by the 
sodar allow production of an image of the structure of the atmos- 
pheric boundary layer. 

It is the purpose of this thesis to calibrate the sodar system 
deployed by the research vessel POLAR QUEEN during MIZEX 84 to 
obtain the  refractive-index structure parameter. The  time-height 
structure of the boundary layer will also be investigated. Efforts 
will also be made to compute the height of the inversion layer and 
to estimate its thickness, the jump of potential temperature, and the 


refractive-index gradient. 
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II. THE EXPERIMENT 


A. LOCATION 
1 Overview 


Ihe recent phase of the Marginal Ice Zone Experiment was 
conducted during the summer of 1984. This experiment, which 
followed that of 1983, was designed to explore the drift, the forma- 
tion, and the evolution of the ice margin through mesoscale 
processes. 

The experiment began with a preparation period consisting 
mainly of satellite imagery analysis. Following this period, seven 
ships, eight aircraft and four helicopters belonging to the ten 
european and north-american countries involved in the experiment 
were deployed in Fram Strait area between Greenland and Svalbard 
Island (Spitsbergen). The first platform was on station at the 
beginning of June and the recovery period started on 17 July 


(Johannessen and Horn, 1984). 


2. Polar Queen Station 


The main goal of the POLAR QUEEN was to folow the 
drifting motion of the ice. Therefore, she was moored to a floe a 
few kilometers inside the ice edge and drifted passively with the 
ice floe while various experiments were deployed to study the 
parameters which described the atmospheric boundary layer. Two 
different drifting phases occurred. The first one ended on 16 June 
when the floe was broken by the propagating swell generated by 
northerly winds. The station was then moved about 60 kilometers to 
the northwest and observations terminated on 17 July. 

This most important part of the experiment was located in 


an area delineated by latitudes 80°N to 81°N and longitudes 001°E 
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to 007°E. We will limit this study to the last week of the experi- 
mentwetrron li July at 21507 to 17 July at 118307 where a most 
interesting data set was obtained. During this period the sodar 
display terminal showed a well formed inversion layer. This 
terminal produced a time-height cross section of the backscattered 
signal using shaded characters whose strength is proportional to the 
signal intensity. This gave a good in situ picture of the ABL 


structure. 
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Fallback Floe, 1 July 1984, Scale 1:7200 





Figure 2.1 The Polar Queen Station. 
The ship was moored on the port side to a large floe of 


approximately 400 by 700 meters as shown in figure 2.1. The 


POLAR QUEEN can be seen in dark, moored along the upper edge 
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of the main floe (named "Fallback floe"). As it Man be seen in 
this picture, the surrounding region is composed of floes of a 
variety of sizes. The floes are all drifting, but not necessarily with 
the same motion. So the environment of the station varied, and this 
can induce variations in the heat and moisture flux depending on 
location. A parcel of air which is at a low height will be under 
the influence of the “underlying ocean surface. A higher parcel can 


be affected by both ice and water surfaces. 


B. WEATHER CONDITIONS 


The weather was characterized generally by stratus or fog 
conditions, although there was less fog than during the 1983 exper- 
iment. Several low pressure systems passed through the region 
during the experiment. For the period of interest, the area was 
characterized by a high pressure system to the east of Svalbard 
Island. Another weaker high pressure center lay over Greenland, 
while a low pressure system was located north of Greenland to the 
northwest of the POLAR QUEEN station. These systems moved 
slowly and erratically. The pressure varied from 1004 millibars to 
1020 milibars, with the minimum values occurring on about 11 July, 
with higher values later. No front crossed the area during that 
week. Winds were variable with velocities ranging from 0 10 193 


m/sec. 


C. MEASUREMENTS 


The two different sources of the data to deal with were the in 


situ measurements and the sodar measurements. 
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1. In Situ Measurements 


In situ measurements were obtained at two different loca- 
tions. One was the ship itself, where surface pressure, air 
temperature, dew point temperature, wind speed and relative direc- 
tion, and downward radiation were measured. 

The second. location was an ice mast located on the ice 
floe at about 50 meters from the edge as shown on figure 2.1 at 
spot "A". It provided air temperature and wind speed at different 


levels, from 5 cm to 5.3 m above the floe surface. 


2 Sodar Measurements 


To avoid noise interference from the ship, the SODAR 
system was installed on the floe at site "B" as is shown in figure 
2.1, about 30 meters from the ice edge. It is expected that the 
acquired measurements were more influenced by surface ice than by 
the surrounding water. 

The sodar system will be described in more detail in the 
next chapter. However, the system gave essentially a backscattered 
Signal which is proportional to the small scale fluctuation of temper- 
ature in a layer encountered by the acoustic signal. These fluctua- 
tions influence the refractive index structure parameter Cn?, which 
is mainly dependent upon the temperature structure parameter CT. 
This parameter is generally larger in the surface layer and in the 
inversion (or interfacial) layer. The sodar gave the values of the 
backscatterred signal for each height at a rate of about one profile 
every minute. 11 wil be possible to deduce the height of the 
inversion layer, and some information regarding its strength may be 
inferred. One goal of this study is to calibrate the sodar by 
computing the value of CT? in the surface layer with only in situ 
measurements and then comparing these values with the sodar signal 


strength to obtain a calibration curve. 
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D. GENERAL DATA REDUCTION. 


Y. Determinati ation of | Surface _ Layer Scales Obtained from 
Surface Layer Similarity — p ج‎ 


The in situ measurements include humidity and wind 
speed and temperature at different levels. From these U*, T* and 
q* may be obtained using: U*=< uw ee EN ٦۷۹۷۷۷2 > وٰ+‎ 0 
g*=<-w'q'>/U*. The bracket represents an average over time. 


These quantities were computed using the bulk method (Businger, 
1973) derived from  Monin-Obukhov similarity theory. A brief 
description of the method and results follows. 

According to the bulk method, the following equations can 


be used to compute U* and T*: 


Us = KUZ ARR) (2210 
T = ka(T-To) / (mZ Z0) > $ (£)) (2 21 
where: 
kK OS von Karman's constant. 
Cl OS ratio of heat transfer to momentum transfer at à - O. 
Ulz? (m/s) wind speed at height ረ. 
ZO, Zo roughness length for velocity and temperature profiles. 


4, )5( and ¢, (5) closed functions of the similarity height parameter 
3 

č is defined as the ratio Z/L where L is the 
Monin-Obukhov length scale (Monin and Obukhov, 1954) 


to: Z/L s Zkg (T* *6.1 10^*Tq*) / T U* (2.3) 


where g = 9.8 m/s? is the acceleration of gravity and q* (gm/kg) 
is the water vapor mixing ratio scaling parameter. 

€ was computed neglecting q* because of the high relative 
humidity and low temperature which produce a low humidity gradient 


and, hence, a low value for q*. Equation 2.3 becomes: 
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ረሺ ን Kg TT Ur (2.4) 


Equation 2.1 can be put in the form Utz? = m f(Z,&) + n, where 
f(Z,5) can be computed using an iteration process to get ደዴ. m = 
U*/k is the slope of the curve and can be measured giving the 
value of U*. n = U*ln(Zo)/k is the intercept and allows the compu- 
tation of Zo. In fact, the drag coefficient was first computed using 


the formula: 
Caves] U/Ut??: (255) 


ros the wind speed at 10 meters. The range of "Cd" 
was found to be between 2.1 and 2.3 10^? m/s. 
T* is compute in the same way as U* and the two results 


produce the following transfer coefficient: 
ሺ. UU z' (T-To) 


where T*/(T-To) = 0.032 is obtained from a plot of T* versus 
(I-To) This leads to: 


CHT 15 1O07 


(Personal communication with P. Guest, 1985). 


2. ኪክ በህ. Considerations for the Calculation of the Heat 
ux 


The heat flux can be computed easily once we know U* 


and T= 


Qo U T 1 )2.6( 

A summary of the most useful data for this study obtained 
from the in situ measurements and averaged over 15 values (about 
122951 


every 3 hours) are presented in Table I . Here the character 


means that there are no reliable data for the corresponding period. 
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IN SITU AND REMOTE MEASUREMENTS 
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A question which arises at this point is: "What was the 
relative influence of the sea and of the ice on the computation of 
U* and T* ?" For the ice mast, there was nearly no influence of 
sea water because the measurements are taken close to the floe 
surface. However for the ship, which was in the water surrounded 
by scattered pieces of ice, the problem may be more complex. The 
surface of the ice can de considered to be nearly at the same 
temperature as the sea water, close to 273°K. It can also be 
assumed that the humidity of the ice surface was nearly equal to 
100 % because of a generally thin skin of melting ice and snow on 
the floes, which more or less resembled large ponds of water in 
this last part of the experiment. Therefore, there will be no 
effort to separate the influences from water and from ice in the 


derivation of U*, T*, q* and Qo. 


3. Calculation of q* 


Assuming than the previous coefficient Ch can be applied 
to the formula which gives q, it is possible to examine the 
assumption that the effects of q* are negligible. An alternative 


way to express q* and T* is the following: 


elo 
^ 


CARO To) (287) 


c'/* (q-qo) (2.8) 


1# 
መሙ 

an 
yl 


where "q" and "qo" are the specific humidity at 10 meters and at 


the surface respectively. "c" is the drag coefficient and is related 
GS Ch by: 

c!/? 2 Ch U(!/U* (2.9) 
If the drag coefficients are the same for q* and T*, q* can be 


expressed by: 
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q* = 1.5 10 UC" qo ٤ (2-10) 


(q - qo) is computed using Clausius-Clapeyron equation (Fleagle 


and Businger, 1980): 


q-qo= exp[(9.4-2353/T)Ln1l0](Rh/100 -1)0.622/P (251) 
Here P is the pressure in  millibars and Rh is the relative 
humidity on the ship. The measurement is assumed to apply to a 


height of 10 meters. The relative humidity at the surface was taken 
to be 100% because of the large ponds of water on the ice floe and 
the melting snow as stated previously. 

The magnitude: for T* was between 1071 °K and 107? °K 
while 6.1 10 *Tq* goes from 107° to 107° g/kg. So this term can, 


in fact, be neglected in equation 2.3 

4. The Roughness Length 
Assuming that Cd is measured close to neutral stability we have: 

Cd =5 K7 IZ 20) (2.12) 
which leads to: 

Zo >۴ 9 +٥ )2.13( 
for Z = 10 meters and Cd = 2.2 10 we get: 

Zo = 2 millimeters. This value is greater than the usual 


quantity for the roughness length, probably because of the 


irregularity of the ice surface. 
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III. SODAR OPERATION AND ANALYSIS 


A. DESCRIPTION OF THE SYSTEM 
1. Basic Principles 


The acoustic sounder operates much like a classical radar 
but emits an acoustic rather than an electromagnetic pulse in a 
narrow beam and receives a signal which comes from scattering of 
the acoustic wave by air parcels at various levels. This backscat- 
tered signal is due to small-scale acoustic refractive index inhomo- 
geneities which are the consequences of local fluctuations of the air 
temperature and velocity. The effects of the humidity fluctuations 
are negligible because they do not affect significantly the sound 
velocity fluctuations, which are directly related to the changes in 
index of refraction (Parry et al., 1972). This will be verified 
later. 

A monostatic sodar system uses an antenna configuration in 
which the transmitters and receivers are  co-located. Thus, the 
received signal is due to direct backscatter. Since the velocity 
variations produce a maximum scattering in the forward direction 
and do not contribute to backscattering, the received intensity will 
not depend on the air velocity fluctuations and wil be mainly a 
function of the local temperature variations. These are described by 
the temperature structure parameter "CT". Since temperature fluc- 
tuations are largest in the surface layer and in the inversion layer 
of an unstable ABL, there wil be maximum backscattered signal 


intensity in these regions. 
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2. Characteristics of the System 


The system emits a 100 ms pulse at a frequency f- 2 kHz 
and then shifts to reception mode for a few seconds. This time is 
variable, depending on the height being probed. The listening 
period is divided into several segments or "bursts" of 4 ms, each 
which will represent a certain distance covered by the signal. To 
have a complete representation of the vertical layer above the 
antenna, the total sounding depth is divided into altitude intervals 
called "gates" which are characterized by a certain number (n) of 
bursts. An average is made for each gate to give the final backs- 
cattered signal. The altitude interval d of the gates is manually 
variable by changing the number n. The value of d is then given 


by the formula: 
d = 4m0 na os (3.1) 


where V is the speed of sound in the atmosphere and is approxi- 


mated by the formula: 


۷۰0۹ 


where T is the absolute temperature of the medium. For T - 273?K 
we get V - 331 m/s, which was used to determine the height 
scale. Consequently, for a gate of 7 bursts, the interval will be 
d- 4.6 meters. There are 200 gates for each profile which can 
cover a height of 25 + (200 4.6) = 3945 meters which will be a 
maximum height of the backscattered profiles. The measurements 


begin at 25 meters because of transmitter ringing and antenna 
near-field effects. 
Another parameter of interest is the wavelength of the 


signal which is given by: 
L = V/E = 16.5 cm 


The most important scattering by small scale fluctuations is done by 
eddies whose size is half the wave length (about 8 cm in our case) 


of the acoustic wave (Parry et al., 1975). 


22 


3. Dependence of the Received Signal on the Layer Stability 


The backscattered signal depends on the small scale inho- 
mogeneities encountered, which are related to the atmospheric 


boundary layer structure and stability. 
8. Unstable Mixed Layer 


Hydrostatic instability occurs when the surface temper- 
ature is larger than the temperature of the layer above. This 
induces upward (positive) buoyancy flux (and heat flux Qo). The 
resultant mixing extends to the interfacial layer. Then the poten- 
tial temperature away from the surface becomes more homogeneous 
and its variations with height are reduced, producing less backscat- 
tering. The highest values of the backscatter signal level will be in 
the surface layer and in the interfacial layer where the vertical 


gradients of temperature are much greater. 
b. Stable Boundary Layer. 


When the surface temperature is smaller than the air 
temperature, there are downward (negative) buoyancy flux and heat 
flux. This tends to establish a less neutral lapse rate above the 
surface and to increase the temperature variance. In this case, the 
small-scale refractive index variations are stronger, and the backs- 


catter signal level will be larger. 


STRUCTURE FUNCTION AND THE TEMPERATURE‏ 71 ے .م 
STRUCTURE PARAMETER‏ 


1. Description 


The refractive index gradient produced by the interaction 
of the acoustic wave with small scale velocity fluctuations is 


expressed by the structure function Dn(r) which is defined by: 


bor) = N(x) - N(x*r)] >. (lec) 
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Here N is the index of refraction at two points in space 
separated by a vectorial distance r. The bracket represents an 
average over all space. Small-scale turbulence theory predicts that 
the structure function value is determined by the structure param- 


eter Cn? and the scalar separation "r" by the formula: 


Dn(r) = Cn?r?/? (333) 


This is valid only if the magnitude of r lies within the inertial 
subrange of the turbulence, so that isotropy is "locally" valid. 
Including both the temperature and water vapor contribu- 
tions, the structure function is related to the temperature structure 
parameter CT?, to the humidity structure parameter Cq? and to the 


temperature-humidity structure parameter Ctq by: 
Cn^-(79 10 *P/T^)^(GmEE II3CIqS*2 2002» (3.4) 


where P is the pressure in milibars and T the absolute temperature 
(Friehe, 1977). 


2. Surface Scaling 


The structure parameters can be related to the measured 
meteorological quantities through Monin-Obukhov surface layer simi- 


larity parameters (Wyngaard, et al., 1971 and Wyngaard, 1973). 


(ep 


T2272 /?£(£) (3.5) 


0 HOL ES (3.6) 


where f(£) is the empirical function found by Wyngaard, et al. 
and updated by Davidson et al. (1978). The quantity "TA" iS a 
constant approximately equal to 0.6 (Fairall, et al. 1980) The 


temperature-humidity structure parameter is given by: 
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Ctq = R T'"Q*Z-? 3A!/?£(Z) (3.7) 


where "R" is the  temperature-humidity correlation parameter 
(approximately 0.8 under unstable conditions). The value of R is 
not well known in the surface layer. Q* is the water vapor density 


scaling parameter (gm/m?) and is related to q* by: 


Where p = 1.2 kg/m? is the air density. 
The Monin-Obukhov length scale was defined by equation 
2.3 and was approximated, neglecting the q* contribution, by equa- 


tion 2.4 , which is repeated here: 


E = Z/L - ZkgT* / T?* (See) 


The function f(§) has different forms for unstable (àí « 0) and 


stable (£ > 0) conditions and is a constant for neutral (Š5 = 0) 
conditions: 7 
ECE ASE ገደ)””/፤ ( < 0: unstable) (999) 
IE) 1.9 ( = 0 : neutral) (3.10) 
HCOS-E1 901 « 2.42771) (£ > O : stable) CD 


The problem of predicting Cn? is therefore reduced to finding 
values of q*, T* and U*. The ratio CT*/Cq* shows the relative 


importance of T* and q*: 


CT? / Cao T*™2/ A "0 (3312) 


This gives: 


DI E rTs2/—0.6 (1.2 q*)^5 1.4 T**/ q** 
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In equation 3.4 there is a factor equal to 3.2 10^? between the 
contribution of the CT? and the Cq? terms. Therefore the contribu- 


tion of Cq? may be significant (more than 5%) if 
3.210፡(1ብ/.ቨ! 6 ጾ 
This leads to: 
l.4 T*/ q**^ = CT Cae o 
Thus the final ےت‎ 15 


T*/ q* < OIE GETS 

For a 10 % contribution, the ratio should be less than 
0.1. This means that q* should be greater than 10 times the value 
of T'* if its influence is significant. It.is known already that we do 
satisfy the necessary condition for this ratio, and the value of Cq* 
can be neglected in equation 3.4. The same result wil be found 
for the contribution of Ctq, even with the imprecision concerning 
DIRAS 

The next step is to use equations 3.9 through 3.11 in 


equation 3.5 to get CT? in the surface layer for every Kind of 


stability: 
CT 5 A OT ILS ኗ < 0 (3-13) 
CTE AT ኣኗ = 0 (3.14) 
CT^s 4 9T]5 7 4 e) ኗ > 0 (3.15) 
These three final equations wil allow to determine the 
value of CT? in the surface layer with only the in situ data. In 


certain circumstances it is possible to measure CT? directly by 
measuring the variance of the temperature at different points in the 


atmosphere, but this was not done for the data for this study. 
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x sodar Data Reduction 
ny The Data Set 


The sodar data, which include sodar characteristics 
and measurements, were recorded on an 8-inch floppy disk and 
subsequently transferred to 9-track magnetic tapes which can be 
used with the NPS IBM 3033 mainframe computer. One tape 
contains the data for the period going from 11 July at 2150Z to 12 
July at 11007 and another tape the remaining data from 13 July at 
02302" to iy July at 1130Z. Unfortunately, there are some gaps in 


the data set due to technical problems during the transfer. 


(1) Time. The date and time are first obtained. 
Some profiles are missing because of computer "interfacing" 
difficulties, and the interval separating two profiles can vary a 
little during each day. This aspect does not change the reliability 
of the measurements corresponding to each time. The only problem 
this generates is in the averaging process, which will be discussed 


in section B.3.b. 


(2) Gate and Altitude. Each profile gives a value 
proportional to backscattered sound intensity at different heights for 
the same time. The tape gives the value of a burst which 15 
always equal to 4 ms in our case. It also provides the number of 
bursts per gate, which allows computation of the altitude for each 
gate. This interval was constant and corresponded to 7 bursts/gate 
for the first tape (11 and 12 July), which gives an interval of 4.6 
meters as shown earlier. For the second tape it varies from 2 to 
6 bursts/gate. This gives an interval between 1.3 and 4 meters 
from equation 3.1. The maximum heights detectable for these cases 
were, respectively, 290 and 820 meters. To simplify the data set 
of the second tape, it was necessary to make a transformation in 
the signal intensity for each profile and each level and to normalize 
the interval everywhere. A 5 bursts/gate interval was chosen 


because it was the value encountered most of the time (and the 
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closest to the other numbers), so that the mathematical 


transformation causes as little imprecision as possible. 


(3) Noise and Signal.. A background noise 
estimate is obtained during each profile and is subtracted from the 
backscattered intensity, after some corrections which depend on gate 
altitude. These final power values become the usable data for each 
level after an averaging process is executed by the system for each 
profile and each gate. 

Additional noise can be caused by aircraft or 
helicopter operations. Because the gain of the system is rather 
high, it is very sensitive to any external noise sources. The air 
flow around the acoustic enclosure can also produce noise. This was 
most likely when the wind speed was higher than 9 m/s. 

When the system detected noise in excess of an 
adjustable threshold, it simply assigned an arbitrary value to the 
backscattered signal intensity of 1 = 32767, whereas usual values 


were between 200 and 600. 


b. The Averaging Processes 


Because of the sometimes important fluctuations of the 
received signal, it was necessary to smooth the response with an 
averaging technique. For this process as well as for the plot of 
the boundary layer, the programs were adapted from those of Lt. 
Mohn (1985). The first step of the analysis was to read the 
profiles by groups of eight and make an average for each gate. 
Then we have a value of the backscatter signal ranging from every 
1 to 10 minutes. This range is variable for a given day and 
increases until 10-15 minutes for the last two days. As discussed 
earlier, the averaging interval is not constant. On a particular day, 
there is not much discrepancy from an arbitrary constant interval. 
However, there are two periods which wil be the object of some 
comments in the time-height cross-section of the boundary layer's 


section. 
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C. REQUIREMENTS FOR SODAR CALIBRATION 


The method to retrieve CT? implies the calibration of the 


system in order to obtain a relation of the form: 
er ED) 


As CT? and I should be proportional, it is expected to find a 
linear function of the first order for f(I). Tom nde CE iS 
computed in the surface layer using equations 3.13 to 3.15 . This 
will be done for the first six levels, whose height goes from 25 m 
to 41.5 m with a separation of 3.3 m given by equation 3.1 , 
when there are 5 bursts per gate. This was not the case for the 
first tape, but after elimination of data which do not fit for this 
computation, only one pair of data remained whose value of CT? 
comes from this tape. 

In the next step a program plotting CT? versus I is used for 
the six different heights. To get good data for each set, the CT’ 
values are averaged for about every half hour. As the values of 
the in situ measurements are already averaged every 10 minutes, 3 
values are averaged, keeping the time of the middle value. Then 
the values of the backscattered signal are read 15 minutes before 
and after this time and averaged to complete the set of two values 
CATA 

The quality of these pairs must be evaluated before generating 


a CT?- f(I) relation, following some criterion: 


* The height of the level considered in the surface layer must 
be less than 1/10 of the height of the inversion layer for the 
equations giving CT? in the surface layer to be still valid. As will 
seen in figure 4.4b, there was an interesting evolution of the 
inversion layer on 13 July. But this does not help very much to 
get favorable sets of data because the inversion layer height is 


always below 250 m from 1100Z to 1700Z. 


e For the same reason, if a second mixed layer forms from the 


ground and so is very low, or if we have a stronger signal at 
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intermediate level, the 6 levels considered wil be probably higher 
than 1/10 this new height. Therefore, these data which correspond 
to a level which is above the surface layer must be rejected. 

For example, on 13 July from 03002 to 1100Z it will be seen in 
figure 4.4b , a second inversion layer whose height varies from the 
ground to 230 meters. The display terminal shows also the same 
kind of phenomenon on 14 July at 03307 when a second inversion 
layer begins to rise up from the surface layer, reaching a height 
of about 150 meters at 0430Z and 200 meters at 05307. Then it is 
completely mixed with the surface layer at 0600Z. Also on 16 july, 
a new layer forms at 0300Z at a height of about 250 meters. Then 
it lowers to 100 meters at 0415Z before joining the surface layer. 

During all these periods, the corresponding sets of values must 
be eliminated when the height of the inversion layer is too low 
(below 250 meters) and cannot fit with the 1/10 rule's requirement 
because the sodar cannot make any measurement below 25 m. 

Another source of problems is the fact that there are some 
gaps in the data set of the in situ measurements made from POLAR 
QUEEN station. This led to missing values for U*. About half of 
these values are missing between 11 July at 2300Z until 14 July at 
230027 Also T* and U* are not known from 16 July 21007 until 
the end of the measurement period. Therefore, it is not possible to 
form pairs of data for the calibration during these intervals of 
time. 

In the backscatter signal data there are also bad data due to 
noise. There are relatively few such cases, and they are generally 
encountered at the first level. For instance, on 14 July from 1315Z 
until 1700Z, about half of the data set at the first level is 
unusable, but there are stil enough values left for the averaging 
process, which is done about every half hour and during which 295 
profiles are read. 

Since CT? and I are obtained in the surface layer for 6 
levels, it is possible to calibrate the system after having deduced 


the function relating CT? and I at these low levels. After that the 
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value of the backscattered signal needs to be used with this rela- 
tion to get the value of CT? everywhere in the atmospheric 


boundary layer. 


D. COMPUTATION OF THE INVERSION HEIGHT 


It is possible to compute the inversion layer height Zi with the 
data set without generating a time-height cross-section. To do that 
it is only assumed that the inversion height occurs when the 
backscattered signal is maximum outside the surface layer. There 
are some difficulties if there are two different inversion layers or if 
the boundary layer is stable and produces a substantial return 


throughout. 


1. The Problem of Determining the Inversion Layer 


The problem of multiple inversions is not too disturbing 
because a second inversion layer is generally weaker and is not 
accounted for in the algorithm. 

The second problem is more troublesome. The strength of 
the signal just above the surface may be stronger than the inter- 
esting signal. In that case, it is necessary to read the data not 
at the first level, but at a higher one, where the strong surface 
layer return is not included. Sometimes this extension from the 
surface reaches higher than the inversion layer itself. "Therefore, 
the first level to be considered must be determined, depending on 
the particular configuration for a given day or part of day. If not, 
the result will indicate either the top of the surface layer or the 
inversion layer, or even a height above the inversion if the first 
level considered is above the inversion. 

The program must be run a few times to determine the 
first appropriate level to be read with regard to a particular 


interval of time. This first level wil be adjusted and changed in 
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order to be always approximately between the surface layer and the 


inversion layer. 


2. The Inversion Layer Height 


Once this difficulty is solved, the height of the inversion 
is computed using two averaging processes. The first one, as 
before, consists of averaging the backscattered signal for each 
level, taking eight profiles over a time of about ten minutes. 

The second process gives the height of the stronger signal 
and is executed by making a running vertical average which spans 
five levels. The stronger value obtained will give the level of the 


inversion layer. 


3. The Inversion Layer Thickness 


Another useful feature which can be computed from the 
sodar data is the inversion layer thickness. This in done in the 
same program that gives the inversion layer height. The averaging 
processes are the same, but instead of searching for the strongest 
vertical signal, we look for the weakest signal below the inversion, 
which should occur where the structure function parameter is the 
smallest. This happens generally at the base of the inversion layer 
Zl. To find the value of the top of the inversion layer Zu, the 
level above Zi is determined where the averaged value of the signal 
is half of its maximum. This factor 1/2 is arbitrarily chosen, but 
since the value of the signal decreases relatively fast above the 
inversion than below, a strong change of the signal wil occur for 
a small height difference. Therefore, the imprecision in choosing 
this factor should not be important. 


The final value Th of the inversion thickness is given by: 


Th = Zu - Zl (3.16) 
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A listing of this program with some comments is given in Appendix 


BA 
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IV. THE RESULTS 


In this chapter, we wil look at the three different kinds of 
results obtained with the methods described previously: the calibra- 
tion of the sodar, the  time-height cross-section history for 
MIZEX-84, and the estimation of the potential temperature and the 


refractive index gradient in the inversion layer. 


A. THE CALIBRATION OF THE SODAR 
1. Analysis of the Results 


Figure 4.1 shows CT? versus I for the six backscatter 


levels closest to the surface. Specifically these are: 


Figure 4.1a : level 1 , 29 meters. 

Figure 4.1b : level 2 , 28.3 meters. 
Figure 4.1c : level 3 , 31.6 meters. 
Figure 4.1d : level 4 , 34.9 meters. 
Figure 4.le : level 5 , 38.2 meters. 
Figure 4.1f : level 6 , 41.5 meters. 


a. Level l 


The first plot seems to be unreliable because all the 
values of I are in the same range, whatever the value of CT?. 
The reason may be related to effects of transducer reverberation on 
the gates which are the closest to the sodar system. This is 
suggested by the fact that a lot of values for I equal to 32767 are 
encountered at this level, an indicator of signal saturation. An 
interesting feature is the asymptotic approach to 1-480 that the 
curve appears to have as CT? increases. Before making any 


comment, the other levels shall be examined. 
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b. evel 2 


This level, shown in figure 4.1b, looks rather like 
the former one. The reverberation still appears to be significant at 
this range, although for small values of CT?, the intensity is lower 
than before. There is evidence of a monotonic functional dependance 
A tor 400212457 (..1 510 *, I=400) Sto ((.1=2 1077. 
[=475). After that point there is, as before, a limitation which 


looks like a saturation of the system at 1=485. 


C. Level 3 


The lower points spread a little more than previously 
and the starting point of the curve seems to be at I-370. 


A few points seem to be out of range: 


. There are five of them in the area delineated by 0<CT*<10%° and 
408 < I < 422. One of these points is taken on 14 July at 1712Z, 
when the inversion layer is about at 320 meters for the half hour 
average. This value is at the limit of ten times the height of level 
3 (31.6 m), and therefore the bulk formula calculation of CT? may 
be tenuous. It may also be due to non-stationarity, in which case 
the fluxes may no longer be constant with height and the similarity 
theory no longer holds. Advection of heat and velocity may play a 
similar role, which invalidates similarity theory by changing the 
balance of the terms in the equations of motion. 

Ons 15 Jule other points span the tune from 09232 
to 10247 and occur at 15287 when the inversion height is between 
310 and 320 meters. As before, these values are just barely in the 
limit accepted for the data set. This situation may explain why CT? 


is smaller than it should be given the value of the intensity. 


e Apart from this group of points, there are also two points 
(CT?z0.6 10^?, I=465) and (CT*=10"*, I=483) which lie well below 
our fitted curve. These two consecutive points occur at 0450Z and 


UC oh 16 July. Before 04157 there was a second and low 
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inversion layer, as stated earlier. After 0415Z, there is no 
apparent reason why there should be any cause to eliminate these 
two values. The disagreement may be due, as before, to non- 
stationarity, advection or unusual ambient noise which lead 10 
stronger values of I compared to real value of CT*. The next 
point, which (CT?-2 10^?, I=479) was obtained at 0626Z is in much 


better agreement with the fitted curve. 


d. Level 4 


At this level, the zero of the curve seems to start at 
about 1=350. There are still two points which appear to be out of 
range, (CT?-0.6 10^, Is459 and  CT*20:96 107? 12478), probably 
for the same reasons suggested in the previous paragraph. 

There are three different domains evident in this plot. 
The first represents the points with low values of CT*, which are 
around a value of I-350. Some of these points have very low values 
of the backscattered signal which occured during the first two days 
of the study period. There is a lot of scatter among the points 
in this area, probably because condition were near neutral, with 
small values of  Ta-Ts (difference between air temperature and 
surface temperature). So the relative effect of instrumental error on 
Ta-Ts may be very large. This scatter may also be due to the fact 
that, even for low values of CT?, the system gives too-high values 
of I because of the background and system noise. 

The second area contains the seven points around 
CT?-0.5 10^? and "5-400. They al correspond to the same period 
from 00002 10 03202 on 16 July, and are described satisfactorily by 
a straight line segment. The third area includes points where CT? 
> 1.5 107? and I > 460. Again, all of these points occur during 
the period from 0615Z to 20407 on 16 July. In this case, the 
points may be described by a straight line from 
(CT*=1.8 107?, Is470) to ስርርር ንበ 1+: ٠٥ Therefore, the 


overall set of points can be described by three different line 
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segments although the first one will not be very reliable because of 


scattering. 
e. Level 5 
Starting at this level, a lot of profiles whose inver- 
sion layer height is less than 380 meters must be rejected. These 


cases are essentially in area one where we got low values of CT? 
in level 4. So little information is lost and most of the interesting 
profiles of level 4 are stil there in that case. 

It can be noted that the slope of the fit looks less 


curved than for previous levels. 


fe Level 6. 


This level gives results similar to level 5, although 
we do not have a well-defined straight line in the second area. But 
the general tendency remains the same and shows good agreement 


with the previous plots. 


2. Calibration of the System 


To obtain a reliable calibration curve, the plots that are 
likely to be free of reverberation effects will be kept. These corre- 
spond to the last three levels (4, 5 and 6). Here the points 
corresponding to the profiles taken at 04507 and 0556Z on 16 July 
are eliminated for the reasons described earlier. The combination of 
these three plots is shown in figure 4.2. 

The three sets of data superimpose quite well but do not 
fit on a straight line as expected if there is proportionality between 
CT^ and I. However, it is useful to plot the same data in loga- 
rithmic coordinates, log(CT?) versus log(I), to investigate whether a 
power relation provides a better fit for the relation between CT* 


and 1. The resulting plot is shown in figure 4.3. The two 


43 


8 قد 1200575 


005 085 09] Ott Oct OOF 082 DS Ort Oct UU 08€ 


ES cp e eb ^ ry E a Y urs LI ال‎ ርጉ) Lj 
/ ,!ተ2 2 aye CLA L EET ROLE LT at E 
. 


LI. 
EL 


59 


ጆር ያ ወርና A ይ ከው ODA E ናቁ መ መድ ሚእ ዓው ሠ 9 «e c Lease xe MES s adu a MEG. e xa e en s ነው ROO O O کرو ا و ےو ہے ے ہے اہ اھ ےا اوھ ود ے ماع‎ AR LS a 





g 19፡57 
G 18407 = o 
p 19.97 =0 

07 






፥ ወ ወ ወ ወ ወ  ፋ=ዓወ =‏ ے۔ 






'u SG*Iy 017 ዩይ”ሄፎ ፡ገካዛከገውዛቫ 9 ፃ S ^y $39^e7 


ome 


لا 


CO 


686. 


رد 


6. 


5 and 


np e Backscatter Intenstt 


CT” 
at Levels 


Figure 4.2 


44 





C) መለ | 

E 128497 

p 18187 - 
0N3937 


O :I4613H 





71 ؛ © 3 0 








5 and 6. 


Log(CT?) versus Log(I) at Levels 4, 


Figure 4.3 


regions with higher  backscattered signal fit rather well on a 


straight line. The least-squares approximation gives the following 


equation: 

Log(CT^) = 8.63 Log(I) - 25.66 (4.1) 
For lower values of I, there is much scatter in the data. This is 
more apparent on this figure. Therefore, the calibration given by 


equation 4.1 should be used carefully for I less than 380. 


3. Error Analysis 


To find the relative error on CT? the value of the vari- 
ance c? is first computed using the approximation: 
go? = 1/N Efloe (CES Clone CP 


where N is the total number of points considered. o is found to 
be equal to 0.01. The relative standard error on log(CT*) is 
given by: 


dilog (CT ))/lor (CREME ۷۷۷٤٤٦ 
whose left hand side is also equal to: 
d(ln(CT?)7ZIn(G'E*) ہہ )ی۔‎ ١ٰ >> +٦ 
These expressions give: 
d(CT?^)/CT* Ese (EEES lO SS) 
Finally, the relative uncertainty on CT? is: 
d(CT4^)/CT-^ "يد‎ 16 +7 


B.  TIME-HEIGHT CROSS-SECTIONS OF THE INVERSION LAYER 


There are several gaps in the time history of the sodar data, 
so the figures showing the boundary layer and its inversion 
obtained from sodar data are decomposed as follows: 

Figure 4.4a : 11 July 2150Z to 12 July 1100Z. 

Figure 4.4b : 13 July 02307 to 24007. 
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Figure 4.4c : 14 July 0000Z to 2400Z. 
Figure 4.4d : 15 July 00007 to 24007. 
Figure 4.4e : 16 July 0000Z to 2400Z. 
Figure 4.4f : 17 July 0000Z to 11302. 


The data are placed in a two dimensional array, one carrying 
the time, the other the altitude (or the gate). The lowest value of 
the backscatter intensity to be plotted is determined by a value 
"BSMIN". This is chosen for each plot such that the isopleths are 
so crowded, and the surface layer and the interfacial layer are 
most apparent. The lower the strength of the inversion layer, the 
lower the value of BSMIN. An incremental value is also chosen so 
that the contour lines are not too tightly spaced. The plot gives 
essentially a contour of the backscattered signal, but the values of 
CT? according to the previous calibration have been added. 

The resulting plots are first examined and their behavior are 
discussed as it relates to environmental conditions. Before doing 
this, it must be recalled that the summer "night" in these latitude 
is much like day. In fact, the downward radiation during the 


night is roughly equal to 2/3 of the value during daytime. 


1. 11 and 12 July 


This period, which includes the night between days 11 and 
12, shows a smooth and continuous decrease in the inversion layer. 
lhis is associated with a more stable layer which characterized the 
night time and a relatively high pressure (1013 mb) on 12 July 
between 12007 and 24007 which caused subsidence and tended to 
lower the inversion altitude. The low value of the chosen BSMIN 
(290) indicates that the strength of the inversion is not very 
great. The background noise for the period is about 230, not very 


far below 290. 
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2s 13 July 


In the second half of the night we see, as before, a 
decrease in the inversion layer. The stability which tends to occur 
contributes to the development of a period of intermittent strong 
backscatter at intermediate levels which rose and mixed with the 
higher one after sunrise. 

Until noon, the wind speed was low and U* was about 10 
cm/s. Then U* increased to 40 cm/s and remained between 30 and 
40 cm/s until 1800Z. This indicates a stronger mechanical mixing, 
which raises the inversion layer through entrainment. Of course this 
increase may also occur as a response to the incoming solar radia- 
tion, which generated more unstable conditions and therefore 
buoyant mixing during the day. 

Dorman (1985) suggests that gravity waves may contribute 
to the kind of movements of the inversion that we can see during 
this day. This, as well as advection, can cause the deepening and 
then the raising of the inversion layer. It is difficult here to have 
a quantitative appreciation of these effects because of the lack of 
supporting environmental data. A higher value of BSMIN (320) 


indicates the strength of the inversion layer for this day. 


3 14 July 


There were no data between 0503Z and 0911Z. This caused 
a jump in the time scale for the plot of the inversion layer. This 
plot has been smoothed by the averaging process. In fact the true 
inversion height goes from 500 meters at 05007 to 390 meters at 
0900Z. 

There was another gap between 20287 and 2203Z, but the 
inversion stayed nearly at the same height during this interval. 
This is why the horizontal scale of figure 4.4c does not show an 
interruption for this short period. 

As before, a stronger backscatter region developed at 


intermediate level from 0330Z until 0600Z. This is not easy to see 


o4 


on the plot because of the lack of data, but in this case it rises 
up to 180 meters then disappears with the early morning. This may 
be caused by large-scale processes such as waves or advection. 


This day is also characterized by a weaker inversion layer. 


4. 15 July 


There are no data for a one hour period between 1620Z 
and 1727Z. The inversion layer went from 265 meters to 245 meters, 
which is a decrease of less than 10$. Therefore there is no need 
for an interruption in the horizontal scale of figure 4.4e, although 
the reader should be aware of this gap period. 

There was not much change during this 24 hour period. 
The signal in the inversion layer has a value greater than 450. 
This and the 13 July data were the two periods where the inver- 


sion layer wil have probably the strongest temperature gradients. 


5. 16 July 


The heat flux was negative (downward) during this entire 
day, and conditions were therefore stable. The temperature on the 
lee mast al 5 5 meters went from 0 22€ at 0000Z to 10 at noon 
and up to 2°C at 2400Z. The temperature just above the ground 
was around 1?C from 1000Z to 1700Z and continued to increase to 
1.7°C at 2400Z. These values of the ambient temperature were 
slightly higher than the other days. 

This increase of the net heat in the boundary layer 
should be responsible for the continuous growth of the inversion 
layer during all the day. It is not easy to see what role the 
subsidence played because we were in a period of high pressure 
(1016 mb), with subsidence expected to lower the inversion layer 
height. Obviously, this was not the case here. The wind speed 
remained between 4 and 6 m/s until noon and then increased and 


stayed between 6 and 9 m/s. This wind caused more mechanical 
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mixing and may also be responsible for the increasing inversion 


layer height. 


6. 17 July 


The wind continued to be relatively strong during this 
day. It stayed between 7 and 10 m/s and is probably the reason 
why the inversion layer continues to behave as it did during the 
day before. There are no in situ data of U and TU for thi- 
period to get the value of the heat flux. The temperature just 
above the ice and at 5.5 meters are very close, so conditions must 
have been near-neutral. This can explain the weak inversion layer 
which tends to vanish and the greater interval in the surface layer 
between the contour lines than on the previous days. lhis was a 
time of formation of an area of stronger backscattered signals at 
intermediate levels around  03-04007 and 0900Z, which decayed 


subsequently. 


C. ESTIMATION OF POTENTIAL TEMPERATURE JUMP AND 
REFRACTIVE-INDEX GRADIENT AT THE INVERSION LAYER 


p The Potential Temperature Jump 


Wyngaard and Lemone (1980) derived an equation which 
shows that it is possible to relate the mean value of the virtual 
temperature structure parameter with the jump of potential virtual 
temperature at the inversion layer, the height of the inversion and 


the virtual temperature scaling parameter: 
<CTV > TOA AS E (4.2) 
where the brackets mean average over the inversion layer and the 


subscript v indicates virtual temperature. lhis is only valid for 


convective cases. However, since there were no strong free 
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convection conditions in the Arctic, we wil consider unstable cases 
for values of Qo greater than 107°? °K m/s. 

Wyngaard and Lemone also show how to estimate the 
potential temperature jump at the inversion, assuming that the 
surface moisture flux and the mean humidity jump can be neglected 
at the top. These assumptions allow use of normal temperature 


instead of virtual temperature. 


Tv = Tu 
Tv = T 


The first assumption is reasonable because q* is very small in this 
experiment, so the vapor flux U*q* will be negligible. The second 
assumption is not so reasonable because the value of the humidity 
jump at the inversion is not known. The presence of fog in the 
area is an indicator of moisture, but it can be expected that the 
cold temperatures (0°C) will produce a generally drier atmosphere 
(due to lower saturation vapor pressure). Radiosonde data will be 
helpful in that case but are not currently available. Therefore, the 
deduction of the inversion layer potential temperature jump should 
be taken as a first guess. With these qualifications, equation 4.2 


can be rewritten: 


AT 2 2 «CT?» Zi?// T* (4.3) 


The only problem now is to compute the mean of CT? in the inver- 
sion layer. Wyngaard and Lemone recommend assuming a peak-to- 
ned Uu NIS CT* Sof 2 The peak value for CT? is easily 
determined because we know Zi and the corresponding value of the 
backscatter signal. Then equation 4.1 is used to get the value of 
CT? at the peak (CTp?), which allows to infer the potential 


temperature jump in the inversion. Equation 4.3 becomes: 


ከ. 21// T 5 
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2. Refractive Index Gradient at the Inversion Layer 


a. Derivation of the Refractive Index Gradient 


The value of the refractive index gradient at the 
inversion layer can be obtained from the equation for refractivity 


(Meeks, 1982): 
N = 77,6 P/T * 3.73 7 (4.5) 


Where P and e are the atmospheric pressure and the water vapor 
pressure respectively in millibars and T is the temperature in “K. 
Assuming that e is negligible which is compatible with the previous 
assumption to derive AT, and taking the derivative with respect to 


Z, equation 4.5 becomes: 


daN 7/7 dZ = m ten alae (4.6) 
Where: 
mes db TP iz 


n = -77.6/T* 
P at the surface was approximately equal to 1005 mb. The ratio 


dP/dZ is given by the hydrostatic approximation: 


8۲۶/٣۰2 = Po RAA (4.7) 
where R is the gas constant. For normal atmospheric conditions R 
= 287 J°K !kg^'. Assuming that T = 273 "we ۲یت‎ ٤7 
mb/km. The inversion height was typically 400 meters, which gives 


P = 960 mb at the inversion and the new value of dP/dZ in the 
inversion layer reduces to -120 mb/km. 
Then the values of the two parameters m and n are: 
m = -34 km”! 
n 


= -0.001 °K ۰۶ 
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Equations 3.16, 4.4 and 4.6 give the value of dN/dZ 


at the inversion: 


dN/dZ = -34 - (CTp? Zi?^?/ T* Th) (4.8) 


Where Zi and Th are in meters and dN/dZ in km'!. 


b. Refractive Index Gradient and Refraction 


Hitney et al. (1961) developed the Integrated 
Refractive Effects Prediction System (IREPS), which is used to 
determine the radar coverage in an area as a function of the 
refractive condition determined by atmospheric parameters as pres- 
sure temperature and humidity. IREPS classification defines the 


following conditions of refraction with respect to dN/dZ: 


Subrefraction: O < dN/dZ 

Normal: 79د‎ > 22 CoO 
Super refraction: - 4 S 
Trapping: 1.11. 


The modified index of refraction M is related to N by 


the formula: 


M = N + 157 Z (423) 


Where Z is the height in km. Therefore we will have trapping 
when dM/dZ «0 or AM «0 for a positive dZ. 


ር. Computed Refractive Conditions 


The values of AT, dN/dZ, AM and some variables used 
to retrieved these values are given in appendix "B". figure 4.5 
show the profile of AM versus time in the inversion layer. The 
trapping conditions happen for the points which are below the hori- 


zontal line (AM=0). Some particular points have been emphasized by 
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a marker. They depicted the largest 10 $ of Qo values. Under 
these conditions, the computation should be more accurate because 
these values are closest to the free convection case. 

The refraction is generally normal until 18007 on 14 
uv  Atter= this period and unti 22007 on 15 July there are much 
stronger values for dN/dZ as well as for AT. But there are some 
important fluctuations in the range of this data and the values 
obtained for AT look sometimes very high for the corresponding 
thickness of the layer which stays between 80 and 150 meters. 
These fluctuations may be due to the different errors and approxi- 
mations introduced in the computation. The more sensitive errors 
are recalled: error on Zl and Zu which determine the thickness, 
error on calibration of the sodar system which is use to get CTp* 
and approximation concerning the mean humidity jump at the inver- 
sion and the 1/2 factor to retrieve the mean value of CT? from its 
peak value in the inversion layer. The major error is probably 
the free convection assumption. The highest values of AT occurred 
very often when Qo is smallest, which means that the situation is 
less unstable and that equation 4.4 is probably not valid, which is 
not surprising. All these features indicate caution is required in 
interpreting these results. However, it is possible to deduce some 
general tendencies during the covered period. There should be more 
refractive conditions during 15 July with some peaks around 02007 
and between 1000Z and 1100Z. This stronger refraction was expected 
in section "B.4" when looking at the time-height cross-section, 
because of the stronger values of the backscattered signal in the 
inversion layer. On the other hand the same expected result is 
not obtained for the 13 July, which shows only two periods (1642 
A 2 aN/dZ smaller than -79 km’. 

The restriction concerning the validity of the computa- 
tion for unstable cases only, does not allow to have any data after 


15 July at 2200Z, because then Qo becomes too small. 
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3. Comparison of some Results with Radiosonde Data 


At the time this thesis was written, there were only two 
radiosonde profiles available. They are showed in figure 4.6 and 
were taken from 13 July at 1100z and 14 July at 1103Z. The left 
profiles represent temperature and dew point with respect to height 


whereas the right profiles depict the wind speed and direction. 


a. Comparison with Radiosonde of 13 July at 1100Z 


The upper and lower levels of the inversion measured 
on figure 4.6 give respectively 370 and 240 meters, which result in 
a thickness of 130 meters. The corresponding AT is equal to 
3.5°K. The program which computes Zu and Zl gave us 366 and 
273 meters at 1057Z and 369 and 204 meters at 1107Z, respectively. 
Therefore, the two values of the upper level are in very good 
agreement. For the lower level and after an extrapolation between 
1057Z and 11077 to find the value at 11007 (about 1/3 of the 
interval), a value of 230 meters is found which is also in good 
agreement with the radiosonde data. 

The values of AT have been obtain with equation 4.4 
by reading directly the values of T* from the in situ data set and 
the values of I and Zi for the other data set. Each of these are 
averaged about every 10 minutes. It is assumed that there is coin- 
cidence when the time between the two data sets is less than 5 
minutes. The value of AT from 1047Z to 1117Z are: 3.42,..2,0, SG 
2.6, the closest being equal to 2.8 at 1057Z. This is fairly good 


agreement with the 3.5°K given by the radiosonde. 


D: Comparison with Radiosonde of 14 July at 1103Z 


Proceeding as before, we look at the set of values surrounding 
1103Z. The first one is at 1100Z (Zu- 382 m, Zl- 336 m) and the 
other one at 1107Z (Zu= 382 m, Zl- 253 m). This gives at 1103Z: 
382 m for the upper level and 300 m for the lower. The radiosonde 
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Figure 4.7 Display Terminal Shaded Plot 
around 1100Z2 on 14 July. 
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measurement gives respectively 510 and 320 meters. The base of the 
inversion seems to be correct, but at the top there is a big 
difference. À reason which can explain this discrepancy may be 
found on the display terminal whose shaded plot is reproduced in 
figure 4.7. after 1100Z, just above the marked dark inversion 
layer, another shaded area can be seen which may be a weak 
prolongation of the inversion layer. This is probably too weak to 
be detected by the program, although the radiosonde responded to 
ita Also, the decay of turbulence above Zi may be an explanation 
of this feature. This less turbulent part is associated with a 
weaker gradient of potential temperature than is the more turbulent 
part around Zi. Of course, this will disturb the computation of AT 
and dN/dZ. In fact, the radiosonde gives a value of 9.1°K while 
the program gives 3.3°K at 1053Z, 1.4°K at 1100Z and 1107Z which 


are out of range of the measured value. 


Two comparisons are not enough to verify the method, 
but as the occurrence of the extended weak inversion layer is not 
common in the data set, we can expect to have a good agreement 
between the computed variables and the radiosonde data in most 


cases if we are in free convection conditions. 
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V. CONCLUSIONS AND RECOMMENDATIONS. 


The analysis of one particular week's data during the last days 
of MIZEX 84 was made for its well-formed inversion layer. This 
analysis allowed finding a power law to calibrate the sodar system 
which was used to study the atmospheric boundary layer. The 
method gives only two periods where the pairs of values CT? and I 
were reliable for this calibration. Both are during 16 July, from 
0000Z to 0300Z and from 0615Z to 2040Z. In the other cases, either 
there were no in situ data, or the inversion layer was too low, or 
we had formation of a second low inversion layer which does not 
allow use of the usual formulas to compute CT? in the surface 
layer. It will be interesting in subsequent analyses to extend the 
method to other days of the experiment in order to see if the cali- 
bration remains in agreement with the relation found in this study. 

This calibration does not seem to hold for values of I lower 
than ۷" ۶۱۹ بب‎ 7٣ However, this is a minor problem 
because the values of CT? in the inversion layer where the backs- 
cattered signal has a value generally higher than 400, were mainly 
used in this study. 

The time-height cross-sections give reliable pictures of the 
boundary layer and especially of the inversion layer. The numerical 
indications provide some complementary quantitative information on 
the value of the temperature structure parameter in the inversion 
layer. 

The calibration relation was then used to deduce values of CT? 
to estimate the potential temperature jump and the refractive index 
gradient in the inversion layer. The error is of order 23 %. Free 
convection formulations indicate that there should be trapping condi- 
tions on 15 July around 0200Z and 10-1100Z. However, since the 


computations are valid only for free convection, these results must 
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be considered with caution because Qo was never very large, about 
107? *K/(m/s) for the trapping condition period. The highest value 
of Qo was around 2 10^? ?K/(m/s) on 13 July between 1400Z and 
1500Z, but the refractive index gradient indicated a normal propa- 
gation condition in the inversion. 

The strength of the backscattered signal in the inversion 
should have lead to a trapping condition, especially on the 13th 
and 15th of July, according to the time-height cross-sections. This 
was confirmed for the 15 July, with the computation of the refrac- 
tive index gradient. It was apparently not true on 13 July, 
although Qo was higher on that day but probably not high enough 
to be a free convection case as required by the theory. | 

The comparison of free convection-determined AT's with some 
radiosonde data was quite good in one case and erroneous in the 
other because of an higher inversion layer top than expected from 
sodar data. To make a better comparative analysis, it wil be 
necessary to obtain more  radiosonde data. This should be a 
productive avenue for future analysis. 

Another source of error of which we were aware too late to 
consider was that the value of I given by the sodar data is in 
fact the summation of the backscattered power and a spreading loss 


correction term AI which is equal to: 
21 - ባነ). ን ۵0۵ 


where t is the time in ms for a pulse to propagate from the 
antenna to the middle of a gate. Therefore Al is a function of alti- 
tude and goes from -44 (for z - 25m) to about -38 when z = 
500m. This represents approximately a 10 % error in the surface 
and in the inversion layer. This should not cause a significant 
change in the overall results, although it should be incorporated 
into further studies using this sodar system. 

This kind of sodar study could be extended to other areas 
with stronger inversions and surface heat fluxes and could be a 
good means, even for an isolated platform, to detect ducting condi- 


tions when a radiosonde is not available. 
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APPENDIX A 


. COMPUTATION OF THE INVERSION LAYER HEIGHT USING SODAR 


O ),:191919/:.1219/91919/9/019191፡(2/091፡፡9/ወ99|2/(/፡፡ወ/9)966,99/0ጃህሇያታያ በ 
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DATAS 


ZIBB JOB (0985 0262) 0 
“MAIN ORG=NPGVM1.0985P 
XEC 0م‎ 


«ፊ 


FORT.SYSIN DD * 


This program reads sodar data from MASS STORAGE 
Date, time backscatter signal. Then it averages 
the signal (INT) among a number "NGATE" of values 
keeping the time of the middle value. 


Then it looks in one profile for the stronger | 
value of backscatter signal (average over "NZ" 
vertical values) which determines the value of the 
inversion layer height "ZI" (if signal is 


greater than a certain value "BMIN"). It also 
determines the lower and upper height of the layer. 

VARIABLES: 

NGATE = No. of profiles to average. 

NVERT = No. of gates to consider. 

IFL = First gate ti read backscatter signal to 
compute "ZI". This allow to eliminate strong 


values encountered in the surface layer. 

B Minimum backscatter signal level to consider. 
NZ Number of vertical value to average. 

BRUIT Average value of the signal over NZ levelS. 

BMIN laximum average value of the signal over . 
BIBI,BMAX : BIBI Determines the maximum backscatter 
signal level (BMAX) in the inversion layer. 


BMOU = Minimum average value of the signal over NZ. 
BMAR I [[ [! ቨ 
ZCI = Corresponding height of BMAR. 

DD = Number of bursts por gate. 

OUT = Out of range variable. 

V - Speed of sound in the boundary layer 

ZI - [nversion sayen height. 

ZG = Height of the lower boundary of inv. layer. 
SIG = Corr BEDNE signal of ZL. 

ZU = Height of the upper boundary of inv. layer. 
SIGU = Corresponding signal of ZU. 

MAI = Level where we begin to look for ZU. 


Thickness of the inversion layer 


THI 2 ZUM- ZC 
INT (NGATE,256) = Value of the signal read in MS 


DIMENSION INT (8,256), BAC(280), ZL1(280), BMAR(280) 


NZ = 5 
NVERT  - 160 
IFL 85 
NGATE = 8 
BBB = 150. 
OUT = -999 
V T a’ 
x Write the label 
WRITE {6 92) 
FORMAT (1X,58('-')/1X,'|DATE TIME| SIGNAL 212 
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13 


ASAS 


107 


MAI 


QOQOOQOO 


1 'SIGNAL ZLOW|SIGNAL ZUP| Dp 58('-')) 


“ok There are two tapes to read in MSS 
tome UNIT 16: 11 12 July, FLAG = 1 
EN UNIT 11: 13 TO 17 July, "ء۶ ئ۲۱۰۲‎ 


LEGE تج‎ E á- e å- x Sg 2] sc II 01i: መመመ መጨ 


ELAG s 
DOZ IR? 
DO = il 
IF. ETAG. E 


¿US 16. d.h co 25983 (INT (I, DIE SG) 
READ X 10, ENDE =999)(INT(I,J),J=1, 256 ) 
CONTINUE Ue 


“= We need to make some adjustements to the first 
"* level to consider to avoid the confusion with 

"* the surface layer. The number below are choosen 
ik after some run and looking at the display 


fous terminal for each period. 


سے سے سے سے سے سس س پس س س  መጨመልመመ on ot nt + ot nb መመመ‏ بے c y ll o. cd‏ جحت س س س د س © © = © = = 


IER ECAG.EQ.1) GO TO 113 

IF(INT(4,4 0 ا‎ 4,9).GT.10) IFL- 40 
٥۳٣۲٣٣۷٠٦٢٣4 4۱۰٦۰٦۸ AND.INT(4,5).GTI.21) IFL- 70 
IF(INT(4,4).EQ.14. AND.INT(4,5).GE.10) IFL- 80 
IF(INT(4,4).GE.1 IFL - 

INTA) GE. 16) IFL 110 

۱" 2) INT 4,9 کے‎ IFL- 135 
IF(INT(4,4).EQ.17.AND.INT(4,5).GT.5) IFL= 120 


0 
000 NCATE 
K= J + IFL - 


በ) ሠጨቐው 


aleale 
MR 


Disregard too big values in the muera. ASS 
"K" means we compute average for the <th" level. 


050211101 K) GT OO ©0112 
808 2 
Ol TOT + )01م‎ EK) ) 
CONTINUE 


EEN NE 0) -GO TO 31 
BA TI - OUT 


GO 
K TOT / FLOAT(N 
CONTI ue = 


69 


AAAA 


(ጋ 


00000 


00000 


70 


NNVERT = NVERT -_4 


IF FLAG. EQ. 5 Der 
d FLAG. NN PET zs 
DO = 1 


ist = Ms Average signal over "NZ" levels. 


DO 50 J = 2 
IET e - d 
IF (BASS EE BN ) GO TO 50 
٢ ብ K).G BIBI) BIBI = BAC(K) 
1 = + 
AVG= AVG + BAC(K) 
CONTINUE 


IF (N.EQ.0) GO TO 55 
BRUIT= AVG/FLOAT(N) 


“i ZI is choosen at the intermediate of the five 
iek levels where the average signal is stronger. 


APR LT.BMIN) GO TO 52 
217 25. «(DD*4./2./1000.*V) * FLOAT((I+1)+(IFL-14)) 
HRN L. BRUIT 
BMAX = BIBI 
MAI = I 


"* Look for the lowest value of the signal 
"* corresponding to a height "ZL1" to determine the 
i lowest boundary "ZL" of the inversion layer. 


(PRUET GT.BMOU) GO TO 55 
BMO BRUI 
BMAR(I) = BRUIT 
ZL) = 22 +(DD*4, /2./1000.*V)*FLOAT ((1+1)+(IFL-14)) 


** Select the value "ZL" of "ZL1" which is below 
= ZT but above 25 meters which is the first level. 


DO 62 II- 1, 280 
LE (211 (111. GE:ZL OR ZI) LE 25 ) CO TO 
SIG = BMAR(II) 


ONTINUE 
IF (BMIN. LE. BBB) GO TO 42 


** Compute the upper height of the inversion layer 


— 7 7 e lm IEEE IE EMG 


TAE E e OUT ) GO TO 42 

BLIN - 7 

DO 75 = MAL, NVERT 
AVG 0. 


NN = 0 
DO 70 JJ = NZ 
KK= “TEL Ne 1 
LF (BAC KK). EQ. OUT 5 GO TO 0 
Eves AVG. + ! BAC(KK) 
CONTINUE 


IF (NN.EQ.0) GO TO 75 


70 


BRUIT = AVG / FLOAT (NN). 
IF (BRUIT. LT; BLIM) GO 


ZU. 25 «(DD* “à b /1000.*V)**FLOAT ((II+1) + (IFL-14)) 


cim CONTINUE 

[s x Estimate the thickness of the inversion layer 
C 

Ta THI = LUT FZL 

© ^* Write date, time, avera max. backscatter 
a "* signal, Zl. SLEPA P 5 GU, ZU, 7 

E 


WRITE ( $ E NTA 1),124,6), BMIN, BMAX,ZI,SIG,ZL, 
„WRITE 18 35 INT 4,1),1=4,6), BMIN, BMAX, ZI,SIG,ZL, 
( im aui tH 7 7 


35 FORMAT (1X, T 1۰۰070 BETO Y" 
F4.0, 1X, F4 0. DX, FA 0; ML U FAO," ', Fa. 0,'MI', 


i 5 54.0: MELO MI 
42 CONTINU 
998 IF لم‎ E 2) GO TO 999 
GO 
999 WRITE. ES 
456 FORMAT 2e 1X,58 ('- 
END. 


/GO. FT11F001 DD DSN-MSS.S0985.MIZE20,MSVGP-PUBAC, 
DCB=(RECFM= ES, LRECL# 1280 BLKSIZE-6400), 
UNIT 3 30V DISP “(OL EEP 

GO.FT16F001 DD DSN=MSS.S0985. [21155 MSVGP-PUBAC, 
DCB=(RECFM=FB LRECL- -1280, BLKSIZE=6400) , 
UNIT TOV DISE KEEP 

GO.FT15F001 DD. DSN= MSS.S0985. Z BT , MSVGP=PUBAC, 

=(RECEM=FB , LRECL=1280, BLKSIZE=6400) , 

UNETS 3330V, DISP- (NEW, CATLG) 
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| مص نل لل لل ل نے ہہ ہم یہ میہ ہہ دہ زیحزہ وہ ہہ ہف وف ہو مم لہہ07ہ٥6ہاےطےلہنہہےہجلہ‏ وو مات نت ١‏ 
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و د س 
| ہے طخ ص حى لمحن مم و اع احم در ون مدر لج ريح رح لم اسم وص ون اه رح رصاح هج ف حم احرص هه ره رص ص بحر وى POOR SAW ST TOW|W OOM COUIM MN DOOR‏ 
FA ANA ON ١‏ لب ص بص پدے اح رم رم بح رس سي TH‏ +ہ لم لہ ہے ہے ہے ہے ہے ہے ہے ہے ہے ہے ہے ANT AA NAMNNNN TH GIO‏ 
t OOO000000000000000000000000000000000000000NNN iii AAOOO :‏ 
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| 10000000000000000000000000000000000000000000000000000000_ 
ل 
NE NN SS RTE EEE‏ = 
١‏ .بت حبص ص جص جص جم جح جم ۰ہ جص بص لہ ہے یہ لح ہے ہے ہے ہے ہے ہے رح ہم ہے رحہے لح (ملمٔہےہےلہم ن ہے ہے ہے لہ ہے ہے رہ نں)ن‪ں )ن اع و ص اجر ہے ہے 
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من لحن‎ © OO NOOO. ميمص احج ام جم رح وص ص بج حي وه ون رص لے لہ ریہ ہے ینہ ےہ ےلم[ہ 6ہ طض ۔۔زن ٦ہی یی نں ۰ی٥0 ح10- یہ۱ زاح‎ 1 
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rtrd mmt ANNI mm AC a a r4 nina r A | | Nate i de 1 | إل رصم رصم يح سم‎ Ir | 
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١‏ الح ص وص صضاحن ن ٠-۲ ۲۹ (۷ ۳0 O0M AT O0 O0 AMNNOUNMOMLVLNONR ቄፎሓቁተጉው፦4ዘር3ሮር ኤ:ኣወቦርሳስምርጋ-ተት LUO MIAO‏ لہ وق لح یے RARO ONO O‏ أسسم وم ل لحن وي حق زا የግኒግ“ር‏ ول 
([ حو جنم وا عمومو فار وو ون MM N‏ على ص تن لح موك ف f‏ ا مو ف 05 وو ف ام حن وو صا ل FALUN‏ رصم ل ل ۷ہ لح لہ لے ہ6 ہف ۳۷ے MO‏ ےی یم ہے ہے |1 ین ٦ن‏ یہ 
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۱ Ec A o INES NINOS ARO NOS SANS عر وک‎ | 
۱ رمجبہےبے مہےہےبہے رہحخرہ ہہ یہ ہےہےم ہہ ہے ہے ہے ہے‎ AO  ےبمیہ ہی ص۸‎ M ہے ہے‎ ۱ 
یت ڑا‎ NN E ooo Ne Oe 
١: ہب ۲ ےم یے بے یلیہ‎ ASE T OFT ትሎ ኡቦኅዓቦግህፕ"1መ”ርት A A AO UN UT ON D LD 00 00 00 DUT TH DO AMMOAOTIDOOONMLWOMADH O MNT OJ CNON i 
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